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ABSTRACT 

Direct collapse black holes (DCBH) have been proposed as a solution to the challenge of as¬ 
sembling supermassive black holes by z > 6 to explain the bright quasars observed at this 
epoch. The formation of a DCBH seed with Mbh ~ 10^“^ Mq requires a pristine atomic- 
cooling halo to be illuminated by an external radiation field that is sufficiently strong to en¬ 
tirely suppress H 2 cooling in the halo. Many previous studies have attempted to constrain 
the critical specific intensity that is likely required to suppress H 2 cooling, denoted as Jcrit- 
However, these studies have typically assumed that the incident external radiation field can 
be modeled with a black-body spectrum. Under this assumption, it is possible to derive a 
unique value for Jcrit that depends only on the temperature of the black-body. In this study 
we consider a more realistic spectral energy distribution (SED) for the external source of ra¬ 
diation that depends entirely on its star formation history and age. The rate of destruction of 
the species responsible for suppressing molecular hydrogen cooling depends on the detailed 
shape of the SED. Therefore the value of Jcrit is tied to the shape of the incident SED of 
the neighbouring galaxy. We fit a parametric form to the rates of destruction of H 2 and H“ 
that permit direct collapse. Owing to this, we find that Jcrit is not a fixed threshold but can 
lie anywhere in the range Jcrit ^ 0.5-10^, depending on the details of the source stellar 
population. 

Key words: quasars: general, supermassive black holes - cosmology: darkages, reionization, 
firststars - galaxies: high-redshift 


1 INTRODUCTION 

Current models of black hole formation and evolution are strongly 
challenged when it comes to explaining the existence of the pop¬ 
ulation of observed bright z > 6 qu asars jMortlock et alj|201 ll: 
IVenemans et al.|[2015l : IWu et alj|2015l) that are believed to be pow¬ 
ered by supermassive black holes (SMBHs) with Mbh ~ 10® Mq. 
Growing these SMBHs from the remnant black holes produced by 
the first generation Population III (Pop III) stars poses a timing 
problem, as this needs to be accomplished within the first Gyr af- 
ter the Big Bang. P op III r emnant mas ses are p redicte d to be low 
IClark et ^ (see e.g. l201lh : IStacv et~^ (see e.g. l2012h . producing 
seed black hole (BH) masses Mbh ~ 10 —100 M0 . To grow these 
seeds by accretion into the observed supermassive black holes re¬ 
quires that they accrete with a low radiative efficiency at close to the 
Eddington rate for « 800 Myr. This naturally requires the persistent 
presence of a gas reservoir to allow steady accretion. Additionally 
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proper treatment of f eedback from this c entral accreting BH must 
al so be accounted for jAlvar ez et alj 2009h . Recent theoretical work 
bv lAlexander & NataraianH2014h has demonstrated that SMBHs at 
z = 6 could indeed form via light BH seeds by super-boosting the 
growth of the seeds in gas-rich nuclear star clusters where the Ed¬ 
dington limit could in principle be circumvented. 


A promising way of avoiding the challenges faced by 
the light seeds invokes a different scenario wherein seed BHs 
form with mass Mbh ~ M 0 from direct collapse of 

primordial ga s as p r oposed by (see EHenstein_&Loe^Jl99^ 
Oh_& Hjumi^ 200 U Bromm_& LoebI 20031 : Koushiappas et alj 


2 OO 4 I : Lodato & Nataraianll2006h . In order to bypass the formation 


of Pop III stars and make direct collapse black hole (DCBH) seeds, 
cooling and fragmentation of gas need to be thwarted in early col¬ 
lapsed dark matter halos. The physical conditions that allow the for¬ 
mation of DCBHs require that no coolants other than atomic hydro¬ 
gen are availa ble in these haloes ( IVolonterill201(]l : lNataraianll201 ll : 
lHaimanll2013h . Such massive DCBH seeds are therefore theorised 
to form in pristine, metal-free atomic-cooling haloes where molec- 
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ular hydrogen (H 2 ) formation can be suppressed. This occurs due 
to the presence of a critical level of radiation in the Lyman-Werner 
(LW) bands of H 2 in the energy range 11.2 — 13.6 eV, produced 
by ext ernal sources in the vicinity of the halo where the DCBH 
forms dPiikstra et alj|200^ : lAgarwal et al.ll20l4 BA14 hereafter). 
Molecular hydrogen can cool the gas down to ~ 100 K, leading 
to fragmentation since the Jeans mass of such a cloud at a den¬ 
sity n = 10® cm“® and temperature T = 100 K is ~ 1000 Mq. 
The suppression of H 2 formation prevents the gas from cooling and 
forming Pop III stars, as the only available coolant in its absence 
is atomic hydrogen that can cool to ~ 8000 K. The Jeans mass 
of such a gas cloud at a densit )0 of n = 10 ® cm“® and temper¬ 
ature T = 8000 K is 10® M 0 . Runaway collapse of this gas 
can lead to the formation of a seed BH in the nuc lear regions that 


can retain up to 90% of the Jeans mass (see e. 

z. Beselman et all 

l2006l:lReean & Haehneltll2009l:lLatif et alil2013l; 

Choi et alj|20l4. 


Most relevant to the final assembly of DCBHs are then the chemical 
pathways that involve the formation and destruction of molecular 
hydrogen in their host haloes. 


The two most critical chemical reactions that control the H 2 
fraction in the collapsing gas are the photodissociation of H 2 and 
photodetachment of 

H 2 -h 7 LW ^ H -f H (1) 

H~ -I- 70.76 —>-11-1-6“ (2) 

where 7 lw and 70.76 represent the photons in the LW band and 
photons with energy greater than 0.76 eV respectively. Destruc¬ 
tion of H“ is critical to the process of DCBH formation as at low 
densities (n < 10 ® cm“®), most of the H 2 formed in the gas is 
produced throug h the following reactions (e.g. iLepp & ShulJl984 
iLepp et ^1200^ 

H -h e H“ + 7 (3) 

H“ -(- H H 2 -h e“ (4) 

Therefore, the value of the specific intensity of the extragalactic 
radiation field in the LW band (defined here as Jlw, the specific 
intensity at 13.6 eV in units of 10“^® erg s“® cm“^ sr“® Hz“®|3 
is the key quantity that controls the efficacy of the DCBH forma¬ 
tion process. This critical value of Jlw that is required to suppress 
H 2 formation sufficiently to allow direct collapse to occur is com¬ 
monly referred to as the critical, threshold value Jcrit. Many au¬ 
thors have previously attempted to determine Jcrit, with the ma¬ 
jority of these studies assuming that the irradiating external source 
can be approximated as a black-body with a surface temperature 
of either 10® K (a T5 spectrum) that is assumed to be representa¬ 
tive of Pop III stars, or lO'^ K ( a T4 spectrum), taken to be rep¬ 
resentative of Population I I stars ( l0mukail|200ll : [^ang et alj 20 li]l : 
IWolcott-Green et alj201 ll . hereafter referred to as KOOl, CS10 and 
WGl 1, respectively). 

However, as we explore in this work, the actual spectral en¬ 
ergy distributions (SEDs) of the first generation of galaxies are not 
particularly well approximated as black-bodies, and recent stud¬ 
ies have shown that using more realistic gal actic SEDs can have a 
significant impact on the value of Jcrit (see ISugimura et alj|^14l 


^ At densities higher than ~ 10® cm“®, the primary channel for H 2 de¬ 
struction is collisional dissociation (KOOl, CSIO) 

^ Note that other definitions of Jlw can also be found in the literature that 
consider the .specific intensity at the mid-point of the LW band or averaged 
over the band. 


Table 1. Summary of the one-zone models in the literature, listing the 
values of Jcrit derived for T4 & T5 spectra and the self-shielding model 
adopted. 


Study 

'^crit 

(T4) 

crit 

{T5) 

Self-shielding 

CSIO 

39 

1.2x10"^ 

DB96 

KS14 

25 

1.4xl0"‘ 

DB96 

KOOll**] 

18.1 

3040 

WGllW 

WGll 


1400 

WGll 

KS14 


1600 

WGll 

SG15a 

18 W 

1630[®1 

WGll 

This work 

18.8 

1736 

WGll 


“ Computed using the KOOl network, but with updated chemical rate co¬ 
efficients; see SG15a, Table 5 
^ WGll: modified version of DB96 
Run 2 in SGlSa 
Run 5 in SGlSa 


Table 2. Summary of the stellar populations considered in this study. 


ISb / CSf 

M* / SFR 

t. 

Z 

IMF 


(MQ)/(M0yr“^) 

(yr) 

(Z0) 


lA/CA 

lOS-io /o.oi - 100 

10 °-° yr 

0.5 

Salpeterl'^l 

IB/CB 

10®-!° / 0.01 - 100 

10°-° yr 

0.02 

Kroupa^^l 


“ slope: 2.35, Mass intervals : 1,100 Mq 
^ slope: 1.3, 2.3, Mass intervals : 0.1, 0.5,100 Mq 


and lAgarwal & Khochfai|[20I^ hereafter KS14 and BA15, respec¬ 
tively). In addition, the value of Jcrit is sensitive to the details of 
the chemical network used to model the gas: changes in the set of 
reactions included in the network or in the rate coefficients adopted 
for them can lead to variations of a factor of a few in Jcrit (see 
iGloveil 2015 - hereafter SG15a - and Glover 2015b). In simula¬ 
tions performed using a T5 spectrum, Jcrit is also highly sensitive 
to the method used to model H 2 self-shielding: studies using the 
original IPraine & Bertoldil i ll9961 DB96 hereafter) self-shielding 
function find values of Jcrit that are an order of magnitude larger 
than those from studies using the modified version of the DB96 
function introduced in WGll dSugimura et alJiToidl) . To highlight 
these differences and dependencies, we summarise the values of 
Jc-rit reported in literature in Table [T] 

In order to accurately capture the role played by LW photons 
in the cooling of gas in these early galaxies, it is therefore necessary 
to use a chemical model that includes all the reactions important for 
determining Jcrit as we do here. We use the best available values 
for the chemical rate coefficients, and properly account for the fact 
that the high-redshift galaxies that produce the photons responsible 
for destroying H 2 and H“ have SEDs that are not simple black- 
bodies. In this paper, we present the results of calculations that im¬ 
prove upon the prior simplifications reported in the literature. Our 
paper is organised as follows. We begin by briefly outlining our 
methodology in Section|2 We present our results in Sectionj^ and 
follow this in Section |4] with a discussion of their implications for 
the efficiency and feasibility of DCBH formation in the early Uni¬ 
verse. 
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2 METHODOLOGY 

In this work we use the one-zone module of the publicly available 
adaptive-mesh-refinement code Enzo to follow the evolution of 
gas to high densities. This one-zone calculation with the up-to- 
date chemical networks for gas cooling (SG15a) tracks the detailed 
gas collapse in primordial haloes. To delineate the potential sites 
of DCBH formation, we model the incident radiation from an ex- 
temal source -first ga laxies in the vicinity- using STARBURST 99 
jLeitherer et al.in~99^ . 


2.1 Basic model 


The one-zone model employed here follows that o f|Q mu ^l l 200 (]|) . 
We set the initial gas temperature to 20, 000 K, and initial gas den¬ 
sity to 10“^ cm“®. The gas cooling threshold at constant density, 
i.e. the lowest temperature the gas can cool to, is set to the H 2 cool¬ 
ing limit 100 K. The specific internal energy is evolved as 


de d 1 . 

dt dt p 

where the pressure is given by 

pkT 

P = -, 

the specific internal energy is 

1 kT 

® = -7-’ 

7a,d - 1 prnn 


(5) 

( 6 ) 

(7) 


7 ad is the true adiabatic index and A is the rate of radiative cooling. 
The cooling rate and associated chemical net work is evolved by 
the Grackl chemistry and cooling solver jBrvan et alj|2014l : 
iKim et al.l |2^ . The Crackle machinery was originally ex¬ 
tracted from the chemistry and cooling network of Enzo, but has 
been updated to make use of the best available chemical and cool¬ 
ing rate data and to ensure that the chemical network contains all 
of the reactions important for determining Jcrit- Details of these 
updates and the reactions included in the chemical network can be 
found in t he Appendix. 

As in lOmuk^ l l200(]tl . the density is evolved as 


dp P 
dt tcoi 


( 8 ) 


where the collapse time-scale, fcoi 


f col - 



(9) 


To account for the effects of H 2 self-shielding, we assume that 
the gas is shielded by an H 2 column density given by 


Nfi2 = ( 10 ) 

where is the H 2 number density and Aj is the Jeans length. 
Based on this value, we then calculate the reduction in the H 2 
photo-dissociation rate using the self-shielding model of WGll. 
We do not account for any reduction in the H 2 photodissociation 
rate caused by the absorp tion of LW photons by the Lyman series 
lines of atomic hydrogen toaiman et al.lll997h . 

To validate our one-zone model, we calculate the values of 
Jcrit corresponding to illumination by external 10^ K or 10® K 
black-body sources. As our chemical model includes all of the key 


® https://grackle.readthedocs.org/en/latest/ 


reactions identified by SGI5a and uses an up-to-date set of chemi¬ 
cal rate coefficients, we would expect to recover very similar values 
for Jcrit. As shown in Table|T] our results agree to within a few per¬ 
cent with those presented in SGI5a. 


2.2 Determining the photo-dissociation (H 2 ) and 
photo-detachment (H^) rates 

Here, we outline our computation of the rate coefficients for H 2 
photo-dissociation and H“ photo-detachment for a more realis¬ 
tic SED. Following BA14 we define the dimensionless parame¬ 
ter, Jlw, as the specific intensity at 13.6 eV, normalised to 10“^^ 
erg/s/cm^/Hz/sr, for a point source at a distance d (in cm). 

dLW = erg/s/cm^/Hz/sr (11) 

TT 47ra^ 

1 

^ erg/s/cm^/Hz/sr ’ 

Where we define a normalised spectrum, L„, as 

L„ = T—— ^ 10”^^ erg/s/cm^/Hz/sr (12) 

E13.6 

which is the ratio of the spectrum Lr, to its value at 13.6 eV, in 
the units of 10~^^ erg/s/cm^/Hz/sr. 

The rate coefficients for H 2 photodissociation and H“ pho¬ 
todetachment can be written as 


kdi = Kdi/Unw (13) 

kde = KdeCnJhW (14) 

where Kde ~ 10“^° and Kdi ~ 10“^^ (see BA15 for 
more details). The dimensionaless parameters a and j3 encapsulate 
the dependence of the rates on the spectral shape of the incident 
radiation field. For a, we have (KO01,BA15) 


‘' 13.6 

1 / 47rL„ 

a = - / —-(Tydiz, 

Kde J nv 
“o-Te 


(15) 


where L„ is as defined previously, (cm^) is the photodetachment 
cross-section (see e.g. Iwishartll 19791 : llohnll 19^^ . vo.jd and v\s,q 
correspond to the frequency limits (Hz) at 0.76 eV and 13.6 eV 
respectively. 

For j3, we have instead the much simpler expression 
(KO01,BA15) 


“la.e 

/ L^di/ 

P = ‘ (16) 

where zzii ,2 & tzia.e denote the frequency limits (Hz) correspond¬ 
ing to 11.2 and 13.6 eV respectively (i.e. the LW band), and 
^vlw = *^ 13 .6 — I'll, 2 . An alternative definition of the parameter 
P = , i.e the ratio of the SED at 12.4 ev and 13.6 eV, has also 

been used in previous works (e.g. KS14, and see I Abel et alJI 19971 
for more details). Our choice is motivated from the fact that the lu¬ 
minosity at discrete energies prod uc ed by different stellar synt hesis 
codes (e.g. iLeitherer et ^Il999l vs iBtuzual & Charloll2003h can 
vary more significantly than the bolometric luminosity in a band. 
We note here, that both definitions are equally valid as long as the 
same convention is adopted throughout a calculation. 
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Figure 1. Density-temperature plots for the collapse of the primordial gas 
in an atomic-cooling halo in the one-zone runs: DC (solid) and non-DC 
(dashed). The solid curve is produced when Eq. [TTlis satisfied by the val¬ 
ues of k(je and k^i passed as inputs to the one-zone model, or in other 
words, when Jlw ~ ■-^crit- The dashed curve is produced when the gas is 
able to cool, i.e. the values of kqi and kdi values are not high enough and 

</lW ^ ^crit- 

To determine the photodetachment and H 2 photodisso¬ 
ciation rate coefficients, we therefore need to specify three num¬ 
bers: a, P, and Jlw- Our procedure for determining these values 
for atomic-cooling halos illuminated by radiation from an external 
stellar population is the same as reported in BA15 and is described 
in more detail in Section[^below. 


3 RESULTS 

We now present the results of the recomputation of the threshold 
intensity required for direct collapse with two key modifications 
compared to earlier work: and updated chemistry, and a more real¬ 
istic assumption for the SED of the irradiating source. 

3.1 A criterion for direct collapse 

We find that there is a region in the kue-kdi parameter space where 
DCBH formation is permitted. To explore this allowed region we 
run our one-zone code for a large number of kde-kdi value sets, and 
identify those for which our one-zone model leads to an isothermal 
collapse at T = 8000 K with densities up to n = 10® cm“® 
(see Fig.|TJ, beyond which collisional dissociation of H 2 takes over 
(KOOl, CSIO). We find that direct collapse occurs when kde and 
kdi satisfy 

kdi = (s“®), (17) 

where 2 = ‘°gio(kie)-B ^ ^ -3.864, B = -4.763, C = 
0.773, and D = -8.154, for kde < 10"®. We plot the 
curve obtained from Eq. [TT] in Fig. We see that for small kde 
( 10"^ s"®), the valuekdi of kdi required for direct collapse re¬ 

mains roughly constant at « 7 x 10"^ s"®. This implies that most 
of the H~ ions that form in the gas are consumed by reaction [4] 
forming H 2 . These H 2 molecules are destroyed via reaction [T] as 
kdi is high enough. Therefore as long as the production of H 2 re¬ 
sulting from the low H" photo-detachment rate can be countered 


Figure 2. The solid curve shows the criterion for direct collapse that we 
have derived from our one-zone runs (Eq. E), while the diagonal lines show 
the combinations of kde and kdi that we get for a T4 spectrum (dotted line) 
and a T5 spectrum (dashed line). For blackbody spectra, changing the value 
of Jlw changes our location along these diagonal lines, with Jlw = Jcrit 
at the point where the lines cross the solid curve. The shaded regions show 
the values of kde and kdi produced by the stellar populations considered in 
Section lT^ computed for an assumed distance of 5 kpc. The region shaded 
with hatched-lines at 45° shows the results from the ISb model and the 
region shaded with hatched-lines at 135° shows the results from the CSf 
model. It is clear' that neither of the commonly-adopted black-body spectra 
provide a good approximation of these shaded regions. 

by kdi, changes in kde have very little influence on the value of kdi 
required for direct collapse. However when kde becomes large, it 
becomes the dominant destruction mechanism for H", thus render¬ 
ing reaction|4]ineffective. Once this occurs, further increases in kde 
strongly suppress H 2 formation, thus allowing for direct collapse at 
smaller values of kdi. 

The diagonal lines in Fig.l^illustrate the relationship between 
kdi and kde for a T4 (dotted line) and T5 spectrum (dashed line) 
respectively. These lines are defined by a = 2000, /3 = 3 for the 
T4 spectrum and a = 0.1, /3 = 0.9 for the T5 spectrum. The 
enormous difference in a for these two spectra demonstrates the 
importance of taking into account the spectral shape of the irra¬ 
diating source while computing the rate constants (BA15, KS15). 
In particular, we see that for the T5 spectrum, the low H" photo¬ 
detachment rate implies that when the criterion for direct collapse 
is satisfied, photo-detachment remains chemically unimportant. On 
the other hand, for the T4 spectrum, H” photo-detachment is very 
important and strongly affects the required value of kdi. It is there¬ 
fore important to understand the effect of changing the irradiating 
source from a black-body to a more realistic SED, as the relative 
abundances of ~ 1 eV photons and LW photons change signifi¬ 
cantly as a function of the mass, star formation rate (SFR) and age 
of the stellar population acting as the source for the radiation. 

3.2 Stellar populations that can lead to DC in their vicinity 

With this understanding of the conditions that lead to direct col¬ 
lapse, we can now ask what sort of realistic stellar populations are 
able to give rise to the appropriate combinations of kde and kdi that 
satisfy the criterion for direct collapse? To explore this, we create 
a grid of stellar populations on the age-mass and age-SFR plane, 
where each point in the grid represents a unique SED. Based on this 
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Figure 3. Stellar populations that allow for DC. Left: results for the ISb model, lA; right: results for the CSf model, CA. The shaded regions in grey indicate 
the stellar' populations that satisfy Eg. llTI for an assumed separation of 5, 12 and 20 kpc (top, middle and bottom panels, respectively) between the atomic 
cooling halo and the in'adiating source. The contours of Jlw at the respective distances are over-plotted in each of the panels. 



Figure 4. Histograms of Jcrit for the ISb model, lA (left panel) and the CSf model, CA (right panel). The histograms are plotted by splitting the stellar' 
populations into ones with f* > 400 Myr and < 400 Myr. The values of Jcrit are obtained by requiring that Eq.[T3be valid in the grey regions in Eig.[^ 
The solid, dotted and dashed lines correspond to the 5, 12 and 20 kpc separations respectively. 


grid of SEDs, we then calculate the corresponding values for kde 
and kde at various assumed distances between the atomic-cooling 
halo (DCBH host) and the stellar population. We outline the steps 
in our approach below. 

Step 1: We use STARBURST99 dLeitherer et alJll99'^ to gener¬ 
ate individual SEDs (not including nebular emission) for an en¬ 
tire set of stellar populations at different stages in their evolution. 
As in BA 15, we consider two extreme cases for the star formation 
history: an instantaneous starburst model (ISb) where all the stars 
form instantly at f = 0, and a continuous star formation model 
(CSf) in which we assume that the star formation rate M* remains 
constant over the lifetime of the stellar population. Eor each stellar 
population, we need to specify the metallicity, along with the to¬ 
tal stellar mass M* (ISb) or the star formation rate (CSf), and 
which denotes the age of the stellar population (time since the on¬ 
set of star formation) in case of ISb (CSf). The values we consider 
are summarised in Table In the ISb models, we generate 1800 
and 3600 SED models for the lA and IB cases respectively. For the 
CSf case, we generate 15550 and 11400 SED models for the CA 


and CB case respectively. This produces a finer grid of SEDs (de¬ 
pendent on galactic properties) as compared to previous work by 
KS14 where they generated 64 and 208 models for their ISb and 
CSf cases respectively. 

We choose a narrow range in metallicity as BA15 and KS14 have 
demonstrated that the values of kde and kdi are far more sensitive 
to the choice of the age and mass (or SFR) than the metallicity of 
the external source. Furthermore, we assume an escape fraction, 
fesc = 1, and an optically thin intergalactic medium. Changing 
these parameters will only reduce the value of Jlw and should not 
qualitatively affect our results. 

Step 2: For a given SED, we then follow the methodology in 
BA 15 to evaluate the rates (Eq. [n&inii by computing: 

• a, / 3 , 

• Jlw is computed using Eq. for assumed physical dis¬ 
tances of 5, 12 and 20 kpc between the putative DCBH host halo 
and the irradiating source. 

The values of kde and kdi that we obtain using this procedure 
are illustrated in Fig.|2 The region shaded with hatched-lines at a 
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45° corresponds to the ISb model, while that shaded with hatched- 
lines at a 135° corresponds to the CSf model. In both cases, we 
assume that d = 5 kpc0 We see that both models produce a simi¬ 
lar range of values for kde since most of the photons contributing to 
kde are produced by long-lived lower-mass stars. On the other hand, 
the UV photons responsible for H 2 photo-dissociation are produced 
primarily by high-mass, short-lived stars, and so the CSf model 
tends to produce much higher values for kdi than the ISb model. 
We find an overlap between the two regions only when the under¬ 
lying SEDs are very similar, which typically occurs only when t, 
is small. For instance, a 10 Myr old stellar population with a mass 
M* = 10^ M 0 that was produced in a single instantaneous burst 
produces a very similar SED to one that forms with a constant SFR 
of M* = 1 Mq yr“^ for the past 10 Myr0 Although most of the 
stellar populations considered here lie below the curve, the reaction 
rates produced by the stellar models that allow for direct collapse 
span three orders of magnitude in both kde and kdi. We also see 
from Fig.l^that the T4 and T5 spectra commonly used when study¬ 
ing direct collapse do not provide a good description of the values 
of kde and kdi produced by realistic stellar populations. An impor¬ 
tant question now arises: does the intersection of the curve and the 
shaded regions in Fig.l^lead to a single well-defined value of Jcrit? 

3.3 There is no unique Jcrit 

In order to explore the implications of Fig. for the value of 
Jcrit, the properties of the stellar populations, i.e. M* or M* & f*, 
that produce photo-dissociation and photo-detachment rates large 
enough to enable direct collapse are show in Fig. The left-panel 
shows the results for the lA SEDs (ISb), which are characterized by 
their age and stellar mass, while the right-panel shows results for 
the CA SEDs (CSf), which are characterized by their age and SFR. 
Qualitatively, lA, IB, CA and CB, lead to the same result, thus we 
only show the lA and CA cases here, and point the reader to the 
Appendix where additional results of IB and CB are shown. The 
three rows (top-bottom) in each panel correspond to separations of 
5, 12 and 20 kpc respectively. The grey shaded region demarcates 
the parameter space for which direct collapse is permitted. Finally, 
the contours correspond to the indicated values of Jlw- We note 
that in every case, the boundary of the shaded region crosses more 
than one contour, indicating that it does not correspond to a single 
fixed value of Jlw- In other words, there is no single value of Jcrit 
in either the ISb or CSf case. We discuss the behavior of Jcrit for 
our two star formation models in more detail below. 

ISb 

At a distance of 5 kpc, young stellar populations with masses as 
low as M* ~ 10® M 0 and f* ~ 5 Myr are able to satisfy Eq.[T 7] 
leading to conditions conducive for DCBH formation. As we in¬ 
crease the separation, the viable ages decreases and the required 
stellar mass increases to compensate for the scaling of the ra¬ 
diation field, increasing by a factor of 5.76 and 16 for d = 12 and 
20 kpc, respectively. 

As previously noted, the shaded region that permits direct 

^ Note that the coiTesponding regions for d = 12 or 20 kpc can be com¬ 
puted easily, since both k^i and kje scale as d~'^. They are omitted from 
the Fig. merely for clarity. 

® Note, however, that the CSf model naturally produces a slightly bluer 
spectrum with increasing age, since there are more high mass stars at any 
given point in the stellar population’s evolution than with the ISb model. 


collapse is not bounded by a single Jlw contour. Young and 
low stellar mass populations with 10^ M* 10®'® M 0 and 

1 f* 50 Myr are represented by Jlw ~ 700. Stellar popula¬ 
tions older than 50 Myr and with M, 10®'® M 0 lie in the region 
where Jlw ~ 700 indicating that in this case, much smaller values 
of Jcrit would suffice. Indeed, the smallest value of Jcrit ~ 0.1 in 
our parameter space exploration is obtained when M* ~ 10^® M 0 , 
ft ~ 300 Myr and d = 5 kpc. 

A range of Jcrit values can be found by imposing Eq. [m on 
the kde-kdi values and systematically lowering the values of Jlw 
in the contours above, till a minimum critical ’Jcrit’ value is found 
for which the equation is still valid. Histograms of these Jcrit val¬ 
ues are shown in Fig. [4] Assuming that DCBHs grow rapidly to 
Mbh ~ 10® M 0 by Eddington accretion with fcdd = 1 from 

2 = 12 (redshift of DCBH seed formation) to 2 = 7, we split 
the stellar populations with an age limit of 400 Myr, i.e. the age 
of the Universe at 2 = 12. We see that for a source located at 
5 kpc, the value of Jcrit ranges from rs. 0.1 — 1000 depending on 
the age and mass of the stellar population. For larger distances, we 
find less variation in Jcrit, although this largely just reflects the fact 
that fewer of the stellar populations in the parameter range that we 
consider produce enough radiation to cause direct collapse. 

CSf 

In the case of the continuous star formation model, there is less 
variation in the value of Jcrit, as can be seen from the right-hand 
panels of Fig. [3] and Fig. |4] Direct collapse is not permitted when 
Jlw 700, which might lead one to infer that Jcrit ~ 700. How¬ 
ever, this is clearly misleading - as there are large regions of param¬ 
eter space (shown clearly in Fig. [3j where Jlw > 700 for which 
direct collapse still remains impossible. Notably, the actual values 
of Jcrit vary by at least an order of magnitude as we change the age 
and the star formation rate of the stellar population illuminating our 
halo. As in the ISb case, we conclude that there is no absolute value 
of Jcrit that we can determine independent of the properties of the 
illuminating source to assess the feasibility of direct collapse. 

It is clear from Fig. [^ & |4] that any attempt to describe this 
behaviour using a single value for Jcrit is inadequate. The critical 
flux depends in a complicated fashion on the mass and age of the 
stellar population and its distance from the atomic-cooling halo. 
Depending on the value of a fixed Jcrit one may therefore either 
dramatically over-estimate or under-estimate the comoving number 
density of atomic-cooling halos that could host DCBHs. 


4 CONCLUSIONS AND DISCUSSION 

The mean background value of Jlw in the early Universe is ex¬ 
pected to rise with the cosmic SFR, reaching a value of 1 at 
2 = 10 (see e.g. KOOl. lCreif & Brommll200^ . Current estimates 
of the critical value, Jcrit, are at least an order of magnitude higher 
than the mean background, which implies th at DC sites will pref- 
erentially be found close to the first galaxies jAgarwal et aljuol^ 
BAM). However, any theoretical model that attempts to explore 
the abundance of such sites is extremely sensitive to the choice of 
Jcrit . For example, an order of magnitude increase in Jcrit can re¬ 
duce th e DCBH number density at 2 = 6 by three orders of mag - 
nitude jDiikstra et alJl20()8l : lLatif et alJl20l3 : lDiikstra et alj|2014) . 
This emphasizes the need to understand and calculate the value of 
Jcrit produced by galaxies composed of Pop II stars in the early 
Universe. 

In this study, we have shown that a Jcrit threshold is not the 
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correct way to identify direct collapse haloes. The quantities that 
determine whether pristine gas in an atomic-cooling halo can un¬ 
dergo direct collapse are the values of the H 2 photo-dissociation 
rate, kdi, and the photo-detachment rate, kde- Traditionally, a 
simple blackbody or power-law spectrum has been used to model 
the radiation from galaxies, resulting in a fixed ratio of these rate 
coefficients. This made it possible to identify direct collapse sites 
using a single value for Jcrit- We have shown here that such an ap¬ 
proach is inadequate when realistic SEDs are taken into account. 
Instead, one must consider the precise age, mass (or SFR), and dis¬ 
tance of the stellar population producing the radiation field in order 
to know whether direct collapse will occur at a given site. Failing 
to account for the detailed shape of the irradiating spectrum could 
lead to a severe underestimation or overestimation of the comoving 
number density of DC sites in the early Universe. 

Comparing our work to the previous study of KS14, we note 
a difference in approach. KS14 assume a fixed value for the rate 
parameter that goes into computing kdi, while we derive the rate 
parameter for kdi directly from the spectrum itself. This inevitably 
leads to a larger variation in Jcrit (especially for the ISb case) than 
that reported by KS14. This could also be due to the larger range in 
stellar masses considered in this study as compared to KS14 (more 
details are provided in sec. l3.2t . 

We expect these results to depend on the choice of IMF. For an 
extremely top-heavy IMF, the results may vary significantly since 
the stellar populations would be predominantly composed of bluer 
stars, thereby pushing the kdi to higher values, i.e. a regime where 
kdi is insignificant. In this work we have neglected the contribution 
of nebular emission lines. These will only serve to add to the total 
radiation budget of the radiation sources. 

The galaxies irradiating the pristine atomic cooling halo are 
also the source of metals that could in principle pollute the DCBH 
site. However in their cosmological hydrodynamical simulation, 
part of the FiBY project, BA14 found that the DCBH sites they 
identified in their volume were metal free and had not been polluted 
in their past by either galactic winds from neighbouring stellar pop¬ 
ulations, or in-situ Pop III star formation, additionally, studies have 
shown that metal mixing appears to be in efficient at these epochs 
dCen & RiguelmelboOSlJsrnith et alj|2015h . 

So far we have discussed the possibility of DC in the vicin¬ 
ity of a single stellar population. However, the presence of multiple 
stellar populations could further aid DC, as the burden of produc¬ 
ing the right photo-dissociation and photo-detachment rates will 
then be shared by multiple stellar populations (BA14). The ideal 
case scenario would be the presence of an old and young stellar 
population nearby, which can give rise to optimal values of kde and 
kdi, respectively. In this case, it makes even less sense to specify 
a single value for Jcrit, since the strength of the radiation flux re¬ 
quired from each population depends not only on the SED of the 
population but also on how much radiation is produced by the other 
nearby source(s). 
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APPENDIX A: CHEMICAL NETWORK 

The chemical reactions used in this study are listed in Table IaTI 
The che mical network is ba sed on the original Enzo network de¬ 
scribed in iBrvan et al.l ll2014h . but we have added two new reactions 
and updated the rate coefficients used for seven others, as described 
below. This modified network includes all of the reactions in the 
minimal reduced model of SG15a that were identified as being cru¬ 
cial for an accurate determination of Jcrit • 


Al New reactions 


We have added two new chemical reactions to the primordial chem¬ 
istry network implemented within Grackle; the collisional ion¬ 
ization of atomic hydrogen by collisions with hydrogen atoms (re¬ 
action 25) and with helium atoms (reaction 26). Reaction 25 was 
included in the study of KOO1, but has been omitted in most subse¬ 
quent studies of the direct collapse model. However, it is important 
to include this process as SG15a show that the additional ioniza¬ 
tion produced by this reaction has a significant effect on the value 
of Jcrit. Reaction 26 was not considered in any studies of the di¬ 
rect collapse model prior to SGI5a but also proves to be impor¬ 
tant, albeit less so than reaction 25. We note that the rate coefficient 
for reaction 25 at the temperatures of interest is highly uncertain 


I Gloveil 2015b 

); we 

iLenzuni et al.l 

199 1|) 


i et al.l 199lh for consistency with the study of KOOl, but 


do not vouch for its accuracy. 


A2 Updated reaction rates 

H~ formation by radiative association (reaction 7) 

The rate coefficient u sed for this reaction within Enzo is taken 
from iHutchinsI (Il976h . but his fit is valid only in the temperature 
range 100 < T < 25 00 K. We have ther efore replaced it with 
the improved fit given in istancil et al.l (Il998h which is valid over a 
much broader range of temperatures. 


H 2 formation by associative detachment of H (reaction 8) 

The rate coefficient for this reaction has rece ntly been mea s ured i n 
the temperature range 10<r<10'‘K bv iKreckel et alJ (l2010h . 
We use the analytical fit that they give to their experimentally- 
determined values. 


H 2 formation by radiative association (reaction 9) 

The rate coefficient previou sly used within Enzo a nd Grackle 
was an a nalytical fit given i n Shapiro & Kand ( Il987h and based on 
data from iRamaker & PeekI lll976h . However, this fit disagrees with 
the data on which it is based by up to 15% at the temperatures of 
interest in the present study. We have therefore repl aced it with the 
improve d analytical fit given in iLatif et all ll20L5h . which agrees 
with the [Ramaker & PeM ( Il976h data to within a few percent. 


Collisional dissociation of H 2 by H (reaction 13) 

The chemical model used in Enzo accounts for the contribution 
made to the rate of this reaction by direct collisional dissociation, 
but not the contribution coming fro m dissociati v e tunn eling, which 
can dominate at low temperatures. iLatif et aT] ( l2Q14h and SG15a 
have shown that omitting this process introduces a factor of two 
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uncertainty into estimates of Jcrit- We have therefore inc l uded it in 
our chemical model, using the rate given in iMartin et alJ lll996h . 


Mutual neutralization of H by tN (reaction 16) 

Enzo uses a rate coeffi cient for this reaction that is taken from 
IPalgamo & Lend ( Il987h . but which disagrees by a factor of a 
few with other de terminations of the rate at low temperatures 
lOlover et alJ i 2006ll. In our mo del, we use instead the rate coeffi¬ 
cient given in Croft et aT] ( Il999h . which agrees well with the results 
of more recent calculations and measurements of the mutual neu¬ 
tralization rate. 


H photo-detachment (reaction 22) 

Our treatment of this reaction is discussed in detail in Sections 
and [3] 

H 2 photo-dissociation (reaction 23) 

For this reaction, we follow the reasoning of KS14 and approximate 
the rate coefficient as 


ItTT-r 0.1/tde. 

The justification for this approach is that in general, the integral 
of the cross-section times the incident spectrum for HJ is around 
a factor of ten lower than for H~, a result which is largely inde¬ 
pende nt of the shape of the incident spectrum fsee lSugimura et al] 
l2014h . In practice, the behaviour of the gas is not parti cularly sen¬ 
sitive to the rate of this reaction (SG15a, Gloveil2015bll . and so this 
approximation does not significantly affect our results. 


APPENDIX B: RESULTS FROM IB AND CB MODELS 

We plot the results (analogous to Fig. [3] and II from the IB and 
CB models in Fie. lBll and lB2l that qualitatively result in the same 
conclusions drawn in sec. 13.31 


Table Al. List of reactions used in this study 


No. 

Reaction 

Reference 

1 

H -1- e“ ^ H+ -1- e“ -1- e“ 

A97 

2 

+ e“ —)■ H + 7 

F92 (Case B) 

3 

He -1- e“ —>■ He+ + e“ -|- e“ 

A97 

4 

He"*" + e“ —)■ He + 7 

BS60 (Case B) 

5 

He+ + e“ ^ He++ -|- e“ -|- e“ 

A97 

6 

He"*""^ + e“ —>■ He"^ + 7 

F92 (Case B) 

T 

H -1- e“ —> + 7 

SLD98 

8* 

H- + H Ha + e- 

KIO 

9* 

H -1- H+ H+ -1- 7 

L15 

10 

H+ -p H ^ Ha -f H+ 

A97 

11 

Ha + H+ ^ H+ + H 

S04 

12 

Ha + e“ —>■ H -1- H + e“ 

DS91 

13* 

Ha + H-s-H + H-l-H 

MSM96 

14 

H“ + e“ —> H + e“ + e“ 

A97 

15 

H-+H-:.H + H-|-e- 

A97 

16* 

H- + H+ H + H 

CDG99 

17 

H- + H+ hJ - 1- e- 

A97 

18 

H^ -1-e- ^ H-1-H 

A97 

19 

H+ -P H- ^ Ha -f H 

DL87 

20 

H + H -p Ha ^ Ha + Ha 

CW83 

21 

H + H-l-H^Ha-l-H 

ABN02 

22* 

H“ + 7 —)• H -1- e“ 

See text 

23* 

H+ -1- 7 -)■ H+ + H 

See text 

24 

Ha + 7 H + H 

DB96, WGl 1, See text 

25+ 

H + H-:.H++e-+H 

LCS91 

26+ 

H + He H+ + e- -1- He 

LCS91 


Reactions marked with an asterisk denote cases where we have updated the 
rate coefficient data compared to the version used in Enzo. The two reac¬ 


tions marked with crosses have been added to Crackle by u s and are not 
includ ed in the original Enzo chemical network described in iBrvan et alj 

boi4i . 

K ey: A97 - lAbel et al.l Jl997h : ABN02 - lAbel et alj J2002h : BS60 
- iBurgess & SeatoiJ Il960t) : CDG99 - ICroft et alj h999t): C W83 - 
ICohen & WestbCTd ll983i)| DB96 
Il987h : DS91 


Draine^_Bertol^ 1^99^; DL87 
Donahue & Shulll 1199 if F92 - 
Latif et alj 

aiT Tl996h : 


iFerland " alj <1993): KIO - iKreckel et alj J2010l): L15 


20151): LCS91 


S04 - ISavin et al 

IWolcot^Gmei^t 


Lenzuniet^ I D99lh ; MSM96 -1^ 


. . . .. . ^^^linet_^ ___ 

J2004bllal): SLD98 - IStancil et alj <19981) : WGl 1 
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Figure Bl. Stellar populations that allow for DC. Left: results for the ISb model, IB; right: results for the CSf model, CB. The shaded regions in grey indicate 
the stellai' populations that satisfy Eq. ll7l for an assumed separation of 5, 12 and 20 kpc (top, middle and bottom panels, respectively) between the atomic 
cooling halo and the in'adiating source. The contours of Jlw the respective distances are over-plotted in each of the panels. 



Figure B2. Histograms of Jcrit for the ISb model, IB (left panel) and the CSf model, CB (right panel). The histograms ai‘e plotted by splitting the stellai* 
populations into ones with t* > 400 Myr and t* < 400 Myr. The values of Jcrit are obtained by requiring that Eg. llTI be valid in the grey regions in Fig. [3] 
The solid, dotted and dashed lines correspond to the 5, 12 and 20 kpc separations respectively. 
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